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ABSTRACT 

To investigate the formation mechanism of relativistic jets in active galactic nuclei and micro-quasars, 
we have developed a new general relativistic magnetohydrodynamic code in Kerr geometry. Here we 
report on the first numerical simulation of jet formation in a rapidly-rotating (a — 0.95) Kerr black hole 
magnetosphere. We study cases in which the Keplerian accretion disk is both co-rotating and counter- 
rotating with respect to the black hole rotation. In the co-rotating disk case, our results are almost the 
same as those in Schwarzschild black hole cases: a gas pressure-driven jet is formed by a shock in the 
disk, and a weaker magnetically-driven jet is also generated outside the gas pressure-driven jet. On the 
other hand, in the counter-rotating disk case, a new powerful magnetically-driven jet is formed inside 
the gas pressure-driven jet. The newly found magnetically-driven jet in the latter case is accelerated by 
a strong magnetic field created by frame dragging in the ergosphere. Through this process, the magnetic 
field extracts the energy of the black hole rotation. 



Subject headings: accretion, accretion disks 
— MHD — relativity 

1. INTRODUCTION 



black hole physics — galaxies: jets — magnetic field — methods: numerical 



Radio jets ejected from radio loud active galactic nu- 
clei (AGNs) sometimes show proper motio n with appar- 



1124-68 are spinning much less rapidly (a — 0.3 — 0.5) 
( Cui, Zhang, fc Chen 199^ ). A similar rapidly rotating 
black hole is also suggested in the AGN of the Seyf ert 1 



cnt velocity exceed ing the apocd of light c ( Pcaraon ct al 



198l | , | Hughos 19& 4 ). The widol}^ accoptod explanation for 



galaxy MCG-6-30 -15 by the X-ray satellite ASCA ( |Iwa- 
sawa et al. 1996). According to recent (nonrelativistic) 



this phenomenon, called supcrluminal motion, is relativis- 
tic jet flow in a direction along the obs erver's line- of-sight 
with a Lorcntz factor greater than 2 (Recs 1966). Such 
relativistic motion is thought to originate from a region 
fcry cloac to the putative auper maaaive black hole whi eh- 



studies of magnetically-driven jets from accretion disks by 
Kudoh & Shibata (1995, 1997a), the terminal velocity of 
the formed jet is comparable to the rotational velocity of 
the disk at the foot of the jet. Further nonrelativistic 



is th ought to power each AGN (Linden Boll 1960 ^ Roe^ 



simulations of jet formation confirm these result s ( Kudc 
fc Shibata 1997b| , |Ouyed, Pudritz, fc Stone 1997| ), except 



1984). On the other hand, the great majority of AGNs 
are radio quiet and do not produce powerful relativistic 
radio jets ( fiees 1984 ). These two classes of active objects 
(radio loud and quiet) are also found in the black hole 
candidates (BHCs) in our own Galaxy. Objects with su- 
pcrluminal jets, such as GRS 1915- H05 and GRO J1655- 
40, belong to the radio l oud class (Mirabel & Rodriguez 
1994, Tingay et al. 1995). Other objects such as Cyg X- 



for the extremely large magnetic field /high jet-power case 
(Meier et al. 1997, Meier 1999) in which very fast jets 
can be produced. The rotation velocity at the innermost 
stable orbit of the Schwarzschild black hole (r = Brg) 
is 0.5c, where rg = 2G'Mbh/c^ is the Schwarzschild ra- 
dius. In addition, it appears that the poloidal magnetic 
field strength in disks around non-rotating black holes may 
be not extremely strong if the magnetic field e nergy den- 



1 and GS 1131 68 arcrolativcly r 



. , J J sity is comparable with that of the radiation (Begelman, 
adio quiet and produce rrr — — — , rrrTrn ft: ^ mi r ^ ■ , 



Uttle or no jet . 



pro 



What causes the difference between the two classes? 
Recent observations of the BHCs in our Galaxy sug- 
gest that the Galactic supcrluminal sources contain very 
rapidly rotating black holes (normalized angular momen- 
tum, a = J/[GAf|jj/c] = 0.9 - 0.95, where G and A/bh 
are the gravitational constant and black hole mass, re- 
spectively), while the black holes in Cyg X-1 and GS 



Blandford, fc Rees 1984 , Roes 1984 ). Therefore, a jet pro- 
duced by MHD acceleration from an accretion disk around 
a non-rotating black hole should be sub-relativistic and 
very weak. In fact, numerical simulations of jet forma- 
tion in a Schwarzschild metric show only sub-relativistic 
jet fiow (Koide, Shibata, fc Kudoh 1999a), except for the 
case when the initial black hole corona is in hydrostatic 



equilibrium rather than free fall (Koide, Shibata, fc Ku 



2 



doh 1998). 

Several mechanisms for relativistic j et formation from 



rotati ng b lack holes have been pr oposed ( Blandford fc Zna- 
jek 1^77, Takahashi et al. 199C). However, up mitil now 
no one has performed a self-consistent numerical simula- 
tion of the dynamic process of jet formation in a rotat- 
ing black hole magnetosphere. To this end, we have de- 
veloped a Kerr general relativistic magnetohydrodynamic 
(KGRMHD) code. In this paper we report briefly on what 
we believe are some of the first calculations of their kind — 
simulation of jet formation in a rotating black hole mag- 
netosphere. 



2. NUMERICAL METHOD 

We use a 3 + 1 formalism of the general relativistic con- 
servation laws of particle number, momentum, and energy 
and Maxwell equations with infinite electric conductivity 
( Thorne, Price, fc Macdonald 1986 ). The Kerr metric, 
which describes the spacetime around a rotating black 
hole, is used in the calculation. When we use Boyer- 
Lindquist coordinates, = ct, = r, — 9, and 



r3 — 



the Kerr metric 



is written as follows, 



^hlicdtf 



i=i 



hf{dx 



i\2 



2h3^3cdtdx-' 



By modifying the lapse function in our Schwarzschild black 
\/l ~ rs/r) to be a 



, and 



hole code (a 

adding some terms of fl^ to the time evolution equations, 
we were able to devel op a KGRMHD code relatively eas ily. 
(See Appendix C in Koide, Shibata, fc Kudoh 1999a| for 
more details on this procedure and the meaning of symbols 
used.) 

We use the Zero Angular Momentum Observer (ZAMO) 
system for the 3-vcctor quantities, such as velocity v, mag- 
netic field B, and so on. For scalars, we use the frame 
comoving with the fiuid flow. The simulation is performed 
in the region 0.75rs < r < 20rs, < 6* < 7r/2 with 210 x 70 
mesh points, assuming axisymmetry with respect to the z- 
axis and mirror symmetry with respect to the plane 2; = 0. 
A free boundary condition is employed ai r — 0.75rs and 
r = 20rs. In the simulations, we use simplified tortoise co- 
ordinates, X = log(r/rH — 1), where th is the radius of the 
black hole horizon. To avo id numerical oscillations, we use 
a simplified TVD method jDavis 1984 iKoide, Nishikawa fc 
Mutel 19961 , [Koide 1997] , [Koide, Shibata, fc Kudoh 1999aD . 
We checked the KGRMHD code by computing Kepler mo- 
tion around a rotating black hole and com paring with an- 
alytic results (Shapiro fc Teukolsky 1983). 



3. RESULTS 

The simulations were performed for two cases in which 
the disk co-rotates and counter-rotates with respect to 
the black hole rotation. Figures la-c illustrate the time 
evolution of the counter-rotating disk case and Fig. Id 
the final state of the co-rotating case. These figures 
show the rest mass density (color), velocity (vectors), and 
magnetic field (solid lines) in < i? = rsin6' < 7rs, 
< z = rcosO < Trg. The black region at the origin 
shows the inside of the black hole horizon. The angular 

^To calculate inner disk rotation cycles, in this paper we always will assume that the inner edge is located at i? = 3rs, regardless of how far 
inward the edge actually has accreted. 



momentum parameter of the black hole is a = 0.95 and 
the radius is th = 0.656rs. The initial state in the sim- 
ulation consists of a hot corona and a cold accretion disk 
around the black hole (Fig. la). In the corona, plasma is 
assumed to be in nearly stationary infall, with the specific 
enthalpy h/pc^ = 1 + Tp/[{r - l)pc^] = 1.3, where p is the 
rest mass density, p is the pressure, and F is specific heat 
ratio and set F = 5/3. Far from the hole, it becomes the 
stationary transonic solution exactly. The accretion disk 
is located at \cot6\ < 0.125, r > td = 3rs and the initial 
velocity of the disk is assumed to be the velocity of a circu- 
lar orbit around the Kerr black hole. The co-rotating disk 
is stable, but the counter-rotating disk is unstable in the 
region R < 4.4rg. Except for the disk rotation direction, 
we use the same initial conditions in both cases. The mass 
density of the disk is 100 times that of the corona at the 
inner edge of the disk. The mass density profile is given 
by that of a hydrostatic equilibrium corona with a scale 
height of Tc ~ 3rs. The disk is in pressure balance with 
the corona, and the magnetic field lines are perpendicular 
to the accretion disk. We use the azimuthal component 
of the vector potential A^jy of the Wald solution to set the 
magnetic field, which prov ides a unifo rm magnetic field far 
from the Kerr black hole (Wald 1974). Here the magnetic 
field strength far from the black hole is 0.3\/poc?, where 
po is the initial corona density at r = 3rs. However, we 
do not use the time component of the vector potential At 
from Wald solution; instead, we use the ideal MHD con- 
dition E-|-vxB = Oto determine the electric field E. 
Here the Alfven velocity and plasma beta value at the 
disk (r = 3.5rg) are va = 0.03c and f3 ~ 3.4, respectively. 

Figure lb shows the state at t = 30ts, where rs is de- 
fined as rs = rs/c. By this time the inner edge of the 
disk has rotated 0.75 cycles, if we assume the edge is at 
R = 3rs.^ Actually, the edge falls toward the black hole 
and rotates faster at i? = 2rs- The rapid infall produces 
a shock at i? = 3. Irs, and the high pressure behind it be- 
gins to produce the jet. This is the same pressure-driven 
jet fo rmation process seen previously i n the Schwarzschild 
case (Koide, Shibata, fc Kudoh 1999a). 

Figure Ic shows the final state of the counter-rotating 
disk case at t — 47rs when the inner edge of the disk 
rotated 1.2 cycles. The accretion disk continues to fall 
rapidly toward the black hole, with the disk plasma en- 
tering the ergosphere and then crossing the horizon, as 
shown by the crowded magnetic field lines near r — 0.75rs. 
The magnetic field lines become radial due to dragging by 
the disk infall near the black hole. The jet is ejected al- 
most along the magnetic field lines. Its maximum total 
and poloidal velocities are the same, v = Vp = 0.44c at 
R — 3.2rs, z — l.Ors. The mass density plot (color) shows 
that the jet consists of two layers. One is an inner, low 
density, fast, magnetically-driven jet and the other is an 
outer, high density, slow, gas pressure-driven jet. The lat- 
ter comes from the disk near the shock at ii = 3. Irs and is, 
therefore, similar to the gas pressure-driven jet of Koide, 
Shibata, fc Kudoh (1998). The former is new and has 
never been seen in the Schwarzschild black hole case. It 
comes from the disk near the ergosphere, and is acceler- 
ated as follows. As there is no stable orbit at i? < 4.4rs, 



3 



the disk falls rapidly into the ergosphere. Inside the static 
limit , the velocity of frame dragging exceeds the speed 
of light (cfis/a > c), causing the disk to rotate in the 
same direction of the black hole rotation (relative to the 
fixed Boyer-Lindquist frame), even though it was initially 
counter-rotating. The rapid, differential frame dragging 
greatly enhances the azimuthal magnetic field, which then 
accelerates the flow upward and pinches it into a powerful 
collimated jet. 

Figure Id shows a snapshot of the co-rotating disk case 
at t = 47rs. The disk stops its infall near R = 3rg due to 
the centrifugal barrier with a shock at r = 3.4rg. The high 
pressure behind the shock causes a gas pressure-driven jet 



with total and poloidal velocities of v 



0.30c at R 



3.4rs, z = 2.4rs. A detailed analysis shows that a weak 
magnetically-driven jet is formed outside the gas pressure- 
driven jet with maximum total and poloidal velocities of 
V ~ 0.42c and Vp = 0.13c, respectively. This two-layered 
shell structure is similar to that of Schwarzschild black 



hole case ( Koide, Shibata, fc Kudoh 1998 ). The centrifu- 
gal barrier makes the disk take much long time to reach 
the ergosphere, which causes the difference between the 
co-rotating and counter-rotating disk cases. 

To more fully illustrate the physics of the jet forma- 
tion mechanism, in figure 2 we show the plasma beta, 
(3 = p/{B'^/2) (color) and the toroidal component of the 
magnetic field, (contour) in the counter-rotating and 
CO- rotating disk cases at t = 47ts. The blue color shows 
the region where magnetic field dominates the gas pres- 
sure; light red — yellow shows where gas pressure is domi- 
nant; and solid contour line shows negative azimuthal mag- 
netic field {B^ < 0), while the broken line the positive 
value {B^ > 0). The toroidal component of the magnetic 
field B^ is negative and its absolute value is very large 
above the black hole in both cases. The field increases to 
more than 10 times the initial magnetic field. This ampli- 
fication is caused by the shear of the plasma flow in the 
Boyer-Lindquist frame du e to the frame dragging effect 
of t he rotating black hole (Yokosawa, Ishizuka, fc Yabukj 
199l |, [Yokosawa 1993| , [Meier 19991 ). Under the simplifying 
assumption that the plasma is at rest in the ZAMO frame, 
the general relativistic Faraday law of induction and ideal 
MHD condition yield, 



dB^ 
dt 



flBr + /2-Se 



(2) 



where /i = c{h^/hi)d{n'i/hz)/dr, and /2 = 
c{h^/h2)d{^j,/h^)/d6. This expression is almost identical 
to that of w-dynamo effect from the field of the terrestrial 
magnetism. Noted that J2 is one order smaller than /i 
when a ~ 1. 

Where does the magnetic field amplification energy 
comes from? It does not come from the gravitational en- 
ergy or thermal energy of the disk, because frame dragging 
effect occurs even when the plasmas of the disk and corona 
are rest and cool. The only other possible energy source is 
the rotation of the black hole itself. Indeed, the increase in 
the azimuthal magnetic field component (eq. (H)) depends 
on the shear of the rotational variable Vl^. We conclude 
that the amplification energy of the magnetic field is sup- 
plied by extraction of the rotational energy of the black 
hole. 



The distribution of the plasma beta (/3) and the az- 
imuthal component of the magnetic field {B^) of the 
counter-rotating and co-rotating disk cases are quite dif- 
ferent. In the co-rotating disk case, they are similar to 
those of the Schwarzschild black hole case. In the counter- 
rotating disk case, the outer part has a positive azimuthal 
component of the magnetic field {B^ > 0), which is caused 
by the counter-rotating disk, and the outer part has the 
high plasma beta. The inner part has a negative azimuthal 
magnetic field {B^ < 0) and low plasma beta. Note that 
the very high plasma beta region (yellow region) is outside 
of the jet; at this point it has almost stopped and eventu- 
ally will fall into the black hole. The negative azimuthal 
magnetic field is caused by the disk around the ergosphere, 
where the disk rotates in the same direction as the black 
hole in the Boyer-Lindquist frame. 

To confirm the jet acceleration mechanism, we esti- 
mate the power from the electromagnetic field, Wem = 
V • (E -I- J X B) and the gas pressure, Wgp = — v • Vp 
along the line, z = l.lrs which crosses the jet foot (Fig. 
3). At t = 47rs, the gas pressure is dominant in the co- 
rotating disk case (Fig. 3b). However, in the counter- 
rotating disk case, the electromagnetic power is dominant 
near the black hole even through the gas pressure power 
is the same as that of the co-rotating disk case (Fig. 3a). 
The magnetically-driven jet in this latter case is acceler- 
ated by the magnetic field anchored to the ergospheric 
disk. The frame dragging effect rapidly rotates the disk 
in the same direction as the black hole rotation, increas- 
ing the azimuthal component of the magnetic field and the 
magnetic tension which, in turn, accelerates the plasma by 
the magnetic pressure and centrifugal force, respectively. 
(A detailed analysis shows that both component of the 
magnetic forces are comparable.) This mechanism of jet 
production, therefore, is a kind of Penrose process that 
uses the magnetic field to extract rotational energy of the 
black hole and eject a collimated outflow from very near 
the horizon. 

4. DISCUSSION 

We have presented general relativistic simulations of 
jet formation from both counter-rotating and co-rotating 
disks in a Kerr black hole magnetosphere. We have found 
that jets are formed in both cases. At the time when 
the simulations were stopped {t = 47rs (53rs), after the 
inner edge of the disk had rotated 1.2 (1.4) cycles in 
the counter-rotating (co-rotating) disk case) the poloidal 
velocities of the jets were v ~ 0.4c (counter-rotating), 
^ 0.3c (co-rotating), both sub-relativistic. In the co- 
rotating disk case, the jet has a two-layered structure: in- 
ner, gas pressure-driven jet and outer, magnetically-driven 
jet. On the other hand, in the counter-rotating case, a 
new magnetically-driven jet has been found inside the gas 
pressure-driven jet. The new jet is accelerated by the mag- 
netic field induced by the frame dragging effect in the ergo- 
sphere. In this case, existence of a magnetically-driven jet 
is not clear outside the gas pressure-driven jet. A longer 
term simulation may show a three-layered structure in- 
cluding the outer, magnetically-driven jet. 

Unfortunately, the counter-rotating (co-rotating) disk 
case could not be continued beyond t = 47rs {t = 53ts) 
because of numerical problems. We have performed one 
other case previously — the infall of a magnetized non- 
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rotati ng; disk into a rapidly-rota ting; black hole ( KoideJ 
Meierl Shibata, fc Kudoh 1999bD. The disk falls toward 



important components of the jet plasma. Simila r detec- 



the bl ick hole more rapidly than the counter-rotating case. 
At later times (after almost two inner disk turns) it deyel- 
oped a relativistic jet with a velocity of w ~ 0.9c (Lorentz 
factor ~ 2). We believe that, if we had been able to per- 
form longer-term simulations here, in at least the counter- 
rotating disk case the magnetically-driven jet also would 
have been accelerated to relativistic velocities (and possi- 
bly the co-rotating case as well). Despite its low speed, 
the magnetically-driven jet in the counter-rotating disk 
case is nevertheless noteworthy because it extracts rota- 
tional energy from the black hole. While t he process is 
simil ar to th e Blandford-Znajek mechanism ( Blandford & 
Znaje lt 1977| ), it appears more closely related to the model 
of Takahashi et al. (1990). In our case, the electromag- 
netic field energy is transformed immediately into kinetic 
energy in the jet. A more detailed analysis and further 
calculations will be reported in our next paper. 

Recently, Wardle et al. (1998) detected circularly polar- 
ized radio emission from the jets of the archetypal quasar 
3C279. They concluded that electron-positron pairs are 



tions in thr ee other radio sources have been made ( Homan 



et al. 1999), which suggests that, in general, extragalactic 



radio jets may be composed mainly of an electron-positron 
pair plasma. The electron-positron plasma is probably 
produced very near the black hole. The mechanism we 
have investigated here, magnetically-driven jet powered by 
extraction of rotational energy of a black hole, is a strong 
candidate for explaining the acceleration of such electron- 
positron jets. 
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Figure Captions 



Fig. 1. Time evolution of jet formation in the counter-rotating disk case and the final state of the co-rotating disk case. 
Color shows the logarithm of the proper mass density; vectors indicate velocity; solid lines show the poloidal magnetic 
field. The black fan-shaped region at the origin shows the horizon of the Kerr black hole (a = 0.95). The dashed line 
near the horizon is the inner boundary of the calculation region. 

At t = and t = 30ts the state of the co-rotating and counter-rotating disk cases are almost identical. However, at 
t = 47rs, while the infall of the disk in the co-rotating disk stops (due to a centrifugal barrier), the unstable orbits 
of the counter-rotating disk plasma continue to spiral rapidly toward the black hole horizon. This difference causes 
the magnetohydrodynamic jet formation mechanisms in the two cases to differ drastically, resulting in a powerful jet 
emanating from deep within the ergosphere. 

Fig. 2. Plasma beta (color) and azimuthal component of the magnetic field i?^ (contour) in the counter-rotating and 

co-rotating disk cases. Solid lines show negative value of B0 and the dashed lines positive values. The two cases differ 
significantly in structure, with the jet in the counter-rotating case originating much closer to the black hole. 

Fig. 3. Power contribution to jet acceleration along the line, z = l.lrs due to the gas pressure (Wgp) and the electro- 
magnetic force (PVem) for both the counter-rotating and co-rotating disk cases. The jet in the counter-rotating disk 
case is accelerated mainly by electromagnetic forces, while that in the co-rotating disk is accelerated mainly by gas 
pressure. Note that, while the power in the gas jet component is comparable in the two cases, the power in the MHD 
jet component is nearly two orders of magnitude greater in the counter-rotating case than the co-rotating case. 
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